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Abstract 

The flagellated bacterium Escherichia coli is increasingly used experimentally as a self-propelled swimmer. To obtain meaningful, 
quantitative results that are comparable between different laboratories, reproducible protocols are needed to control, ‘tune’ and 
monitor the swimming behaviour of these motile cells. We critically review the knowledge needed to do so, explain methods for 
characterising the colloidal and motile properties of E. coli cells, and propose a protocol for keeping them swimming at constant 
speed at finite bulk concentrations. In the process of establishing this protocol, we use motility as a high-throughput probe of 
aspects of cellular physiology via the coupling between swimming speed and the proton motive force. 
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Some time ago, our lab wanted to culture motile bacteria as 
‘model active colloids’. We obtained a strain of Escherichia 
coli with the full complement of motility genes and a cultur¬ 
ing protocol from a local microbiologist. For some time, we 
thought we were experimenting with motile E. coli, until one 
day we checked in the microscope. Few, if any, of the cells were 
swimming! So we set out to learn how to modify the standard 
protocol to optimise motility by collating literature, talking to 
other researchers and trial and error; we also implemented dif¬ 
ferential dynamic microscopy (DDM) to quantify motility. 

This article reviews what we have learnt. Some of the ma¬ 
terial is previously known, but seldom critically discussed in 
one place. We have explained some basic bacterial bioener¬ 
getics and genetics, because physical scientists can use E. coli 
and collaborate with biologists more effectively if these top¬ 
ics are understood. Much of the materials is new, arising from 
using DDM to quantify motility. While we aim primarily at re¬ 
searchers working on active colloids Q, this article should also 
be useful to others studying motility biophysics |2]| . 

From the outset, we refer to various culture media (BMB, 
TB, LB) and protocols, and freely use terminology related to 
molecular biology (plasmid, gene names, etc.) and measure¬ 
ment techniques (OD, DDM, etc.). Readers should refer to Sec¬ 
tioned! matters of cell culture, Sections[3]and[4]for methodol¬ 
ogy, and Section 10 and Appendix C|for biolog ical jargon. We 
also provide a table of symbols in | Appendix F| 


tive colloids are out of thermal equilibrium even without exter¬ 
nal driving. Unusual phenomena displayed by such systems, 
such as ‘negative viscosity increment’ 0 and ‘rectification’ 
l6l , pose a ‘grand challenge’ to statistical mechanics ]7][8]. The 
self assembly of active colloids, possibly mixed with passive 
particles, may provide routes to new ‘smart materials’ El- 

Progress in this new area will be greatly facilitated by exper¬ 
imental data from ‘model systems’ that can most ‘cleanly’ con¬ 
front theory and simulations. Historically, well-characterised 
model passive colloids have enabled progress at critical points, 
from sedimentation equilibrium maim to hard-sphere crys¬ 
tallisation m and sticky-sphere glass transitions f!3l . The 
minimal requirements for a ‘model passive colloid’ include 
reproducible synthesis/preparation, known particle size/shape 
distributions, quantifiable and ‘tuneable’ interparticle interac¬ 
tions and accurate particle volume fractions M El US El. 
For model active colloids, we also need propulsion mechanisms 
that are understood and ‘tuneable’, and knowledge of motility- 
specific interparticle interactions, e.g., via hydrodynamics. By 
these criteria, there is as yet no ideal model active colloid. 

Escherichia coli is widely used in active colloids research. 
Its self propulsion is understood in essence lfl8l [1911201 , while 
propulsion mechanisms in many synthetic swimmers are still 
debated Enm. However, for E. coli to become a fully- 
fledged model, reproducibility, control and characterisation 
need to be improved. We address these issues in this article. 


1. E. coli as a model active colloid 

An ‘active colloid’ III is a suspension of w 5 nm-5 pm par¬ 
ticles that consume ‘fuel’ to propel themselves, such as motile 
bacteria or various synthetic ‘colloidal swimmers’ 011 . Ac- 
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2. E. coli as a bacterium 

Escherichia coli l23l is the best understood living organism 
on earth today. It is often the model of choice for understanding 
molecular biological processes: ‘Tout ce qui est vrai pour le 
Colibacille est vrai pour T elephant.’ l2l 
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Table 1: Average length and standard deviation of 30 cell bodies of E. coli 
AB1157 grown in LB/TB at 30 °C/37 °C to late-exponential (OD = 0.5) or sta¬ 
tionary (16 hours growth) phase. Note that these cells have not been washed in 
BMB. Washing, however, makes little difference. 
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Figure 1: Visualisation of both body and flagella bundle of swimming E. coli 
from a real-time movie. Flagella are fluorescently labelled using Alexa546 
and the body with a plasmid expressing GFP (ADI pHC60, TableUs¬ 
ing epi-fluorescence microscopy these can be imaged into separate channels 
at high enough frame rates to avoid motion blurr. Image part of a 100 fps 
movie acquired using an Orca Flash 4.0 camera on a Nikon Ti microscope with 
PA60x/1.4 OIL objective using a Cairn Optosplit II. 


The cytoplasm of E. coli is enclosed by two lipid bilayer 
membranes. Between them is a peptidoglycan ‘cell wall’ 
(or ‘periplasmic space’), a network of poly-sugars (‘glycans’) 
linked by short peptides. Within the species known as E. coli 
there are many genetic variants, or ‘strains’ (so that ‘bacterial 
species’ is a problematic concept (25l ). The genome of strain 
MG 1655 was first to be sequenced (26l : today, the genome of 
2085 strains are available in GenBank (Z7l . although > 30% of 
the genes remain of unknown function. The best single source 
of information on E. coli (and its close relative Salmonella) is 
a two-volume ‘bible’ If28l now updated digitally 1291 . 

Each cell carries multiple helical flagella lf30ll3T) powered by 
rotary motors E [T8]| embedded in the membranes and cell wall 
l32ll . When all flagella rotate counterclockwise (CCW) (viewed 
from behind), they bundle to propel the cell forward. When 
a motor reverses to clockwise (CW), its flagellum unbundles 
and the cell tumbles. Wild-type (WT) cells run (« 1 s), tumble 
(~ 0.1 s), and run, the latter in a direction that is more or less un¬ 
correlated with the nin before the tumble (33]. When successive 
sensings of the environment return an increasing concentration 
of an ‘attractant’, intracellular signals reduce the rate of CCW 
—> CW —» CCW switching (i.e. the tumbling rate). This biased 
random walk up an attractant gradient constitutes chemotaxis. 

A number of E. coli strains used for active colloids and re¬ 
lated research are listed in Table B.2 Appendix B Many of 
these were derived from K-12, a ‘debilitated’ laboratory strain 
(341135 1 that does not normally infect humans. It is deemed safe 
for laboratory use and classified as a category 1 biohazard (36]. 
Mutants in which one or more key chemotaxis genes have been 
deleted, e.g. AcheY , lose the ability to tumble. These ‘smooth 
swimmers’ simply swim forward, with a persistence time set 
by their inverse rotational diffusivity D r = ImjLcr /k B T « 5 s, 
where we have modelled a cell body with its flagella bundle as 
an ellipsoid of major and minor axes L « 10 pm and cr « 1 pm. 


For mostly unknown reasons M BD, bacterial shapes are 
highly conserved 0 shape mutations tend to be lethal. In this 
section, we review the single-particle ‘colloidal’ properties of 
E. coli and how to estimate cell concentrations. 

3.1. Shape and size 

3.1.1. Microscopy 

The shape and size of E. coli cell bodies can be obtained us¬ 
ing microscopy and scattering (compared in |[43l ). Measuring 
cell bodies using phase-contrast microscopy of strain AB1157 
harvested using our standard protocol and washed into BMB 
gave (l) = 2.4 ± 0.6pm, where the uncertainty is a polydis- 
persity. At mid-cell, (cr) = 0.86 ± 0.07pm, where the uncer¬ 
tainty represents measurement errors. Importantly, these ‘vital 
statistics’ depend on strain, growth conditions and the time of 
harvest. Table [T] Time of harvest is important because bac¬ 
teria grow in well defined stages: from rapid division in an 
exponential phase (cell number density n grows exponentially 
with time), through a stationary phase (n = constant), to a death 
phase (n decreases with time). Our data show that stationary- 
phase cells tend to be less polydisperse in length than mid- 
exponential cells, as found previously (44]. There may also 
be some dependence of cr on / ga¬ 
it is often desirable to visualise the flagella bundle, because it 
renders the cell highly anisotropic and introduces strong steric 
effects (46)03. If the cell body and flagella are labelled with 
the same dye 1 481149 1. fluorescence from the former will dom¬ 
inate. To visualise the flagella alone, the fliC gen (0c an be mu¬ 
tated to facilitate binding to Alexa Fluor C5 maleimide (e.g., 
Alexa546) (50|. If cells carry GFP-plasmid pHC60 (Table[BT2]i, 
then the cell body and flagella can visualised simultaneously. 
Fig. □ Note that such mutants should always be checked for 
potential changes to their motility and other phenotypes. 


3. E. coli as a colloid 

The cell body of an E. coli bacterium is a net negatively 
charged spherocylinder (diameter cr, pole-to-pole length Of! 


1 For physicochemical 'vital statistics’ of E. coli, see I37H381I391 


3.1.2. Scattering 

Cell shape and size can also be studied by static light scatter¬ 
ing (SLS) either averaged over cell populations (5T1 l52ll53ll54l 


2 It is possible to turn bacteria into silica colloids EU, opening up a novel 
route to synthetic colloids with many shapes (40). 

' This encodes flagellin, the protein units used to build up flagella filament. 
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I55ll56l or individually in flow cytometry l57l . Given informa¬ 
tion on the refractive index (58), models can fitted to scattering 
data. Typically, an E. coli cell with cytoplasm surrounded by 
a cell wall is modelled as a core-shell ellipsoid |51, 59. 60j, 
giving interference fringes that survive polydispersity in / be¬ 
cause there is very little variation in either cr or periplasm thick¬ 
ness t52\Ml How optical inhomogeneities due to the chromo¬ 
some affect dynamic light scattering, which can also give in¬ 
formation on size and shape, has been studied l62l . Polarised 
scattering gives extra information Il63ll64l . 


3.2. Sedimentation 

The (mass) density of E. coli cells, 1.08 <P<1 .16gcm 3 , 
depends on media and conditions |[65l [66. |671. The density dif¬ 
ference with the solvent, Ap < 0.1 gem -3 , causes cells to sedi¬ 
ment. The distribution of particles in a suspension in sedimen¬ 
tation equilibrium defines the colloidal length scale ED. In this 
dynamic equilibrium, concentration-driven diffusive flux from 
the bottom is balanced by a sedimentation flux from the top, so 
that when the suspension is dilute enough to neglect interparti¬ 
cle interactions, the number density profile with height is 

n(z) = n(Q)e~ zlt \ (1) 


where the gravitational height t g is given by 


Do 

v s ' 


( 2 ) 


Here, Z) 0 is the free-particle diffusivity, and v s is the single¬ 
particle sedimentation speed. A colloid is suspendable against 
gravity, so that l g has to be larger than its characteristic size. 

The friction coefficient for motion parallel to the long axis of 
a prolate ellipsoid with major and minor axes / and <j is l68l 


$1 = 



(3) 


where r/ is the solvent viscosity. For E. coli in typical motility 
experiments, l/cr « 2 to 3, so that A « 2mjl, from which 


v 


Apgo- 2 
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(4) 


Translational diffusion along the long axis is governed by 


k B T 

k ’ 


(5) 


so that the gravitational height of a non-motile cell is 


4 c °ii) ~ (±I<T ~ 4 pm * 21. (6) 

ncr~IApg 

Thus, the cell body of non-motile E. coli is just colloidal. 
Motile E. coli should have significantly higher ( li69l : and 
motility coupled with growth can display rich phenomena ITOl . 



Figure 2: Fraction of bacteria sticking to the lower capillary wall in a « 4 cell 
width layer at different concentrations of TWEEN 20 surfactant in BMB with 
0.72mM glucose for WT AB1157 at « 6 x 10 7 cells/ml. Points are averages 
over 0-30 minutes, 30-60 minutes and 60-90 minutes after loading. Arrows: 
adhering fractions with no surfactant. 


3.3. Interactions: stability and surface adhesion 

The surface properties of E. coli ED EE] ED are strain de¬ 
pendent |74| . Interpreting electrophoretic mobilities (p) is non¬ 
trivial because surface macromolecules Gama ns and the 
periplasm are ion-permeable. Even with a ‘soft’ electrokinetic 
theory allowing for this 1751 1761 . relating p to zeta potentials 
and surface structures is fraught. 

Nevertheless, E. coli definitely carries a net negative charge 
under physiological conditions (-3 < p < -2pm cm V -1 s -1 ) 
largely due to ionised carboxylate and phosphate groups lf74l 
ED- This charge is normally sufficient to confer colloidal stabil¬ 
ity in growth media or in BMbQ H owever, a strain that has be¬ 
haved stably for many months may suddenly start aggregating 
due to ‘phase variation’ in the expression of the surface protein 
Ag43 (antigen 43) 1771 . Each generation, there is an 0(10 3 ) 
probability that a cell switches between ‘on’ and ‘off’ states of 
Ag43 expression. ‘On’ renders cells ‘sticky’ and leads to ‘au¬ 
toaggregation’. The expression state is heritable and reversible. 

For active colloids work using E. coli close to surfaces lf78l 
EOTim it is important to minimise adhesion, which is far 
from understood. Various non-ionic surfactants prevent E. coli 
from sticking to untreated glass 148] 1ST) . Figure [2] quantifies 
this effect for a popular surfactant used for this purpose. 

We imaged AB1157 cells in BMB within an optical slice 
of ~ 4 cell width next to the bottom surface of an untreated 
glass capillary using a Mikrotron MCI362 camera and a Nikon 
Ti microscope with 60x phase contrast objective (PF60x/0.7). 
Automated software for identifying Il82l and tracking rods was 
used to analyse movies, classifying cells into ‘free’ (diffus¬ 
ing or swimming) and ‘adhering’ (stuck). The average frac¬ 
tion of adhering cells, extracted from many movies, dropped 
abruptly to > 0 at a TWEEN 20 concentration of ~ 10 2 wt.%. 
Note that concentrations up to 0.2% have been used previously 
ED- Preliminary DDM data suggest that commercial TWEEN 
may contain impurities rapidly metabolisable by E. coli, so that 


4 The screening length in BMB is < 1 nm. 
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high enough concentrations of TWEEN may significantly alter 
swimming behaviour. 

3.4. Cell concentration 

The definition of cell concentration and its determination are 
non-trivial. The cell body volume fraction, <fi, gives a better in¬ 
tuitive feel for the degree of ‘crowding’ than the number density 
of cells, n. The effective volume fraction, which takes into ac¬ 
count the cell consists (body length l ) and flagella bundle with 
total length L < 10 pm is 0 eff ~ (L/l) 3 cf> > 10 2 (f) if the cells are 
randomly oriented, so that they reach ‘overlap’ at <l> < 0.01. 

The standard way to obtain n is from the spectrophotometric 
optical density (OD) of cell suspensions, which is proportional 
to n for n < 10 9 cells/ml (corresponding to <p eS ~ 0.2)J^] Cali¬ 
bration is by ‘viable plate counting’ (Section [53j i, and is found 
to depend on growth stage and other conditions. 
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Figure 3: (a) Fitting DDM data g(q,r), Eq. [7] and [8] using the ISF given by 
Eq.^at one q for WT AB1157 at 5 x 10 8 cells/ml. Main plot: data for bacteria 
swimming fast enough for ballistic and diffusive motions to show up clearly as 
two processes, separated by a point of inflection at r « 0.4 s. Inset: the bacteria 
have slowed down to the extent that swimming and diffusion are no longer well 
separated, (b) Fitted speed v (main plot) and non-motile fraction y 3 (inset) for 
the data shown in the main part of (a) (black) and in the inset of (a) (red). 



4. Characterising E. coli motility 

We turn next to characterising E. coli as self-propelled par¬ 
ticles. The average swimming speed, v, and the fraction, / 3 , of 
non-swimmers that only diffuserlcan be measured using a high- 
throughput technique, differential dynamic microscopy (DDM) 
Il84l l85l l86l . Application to swimming algae f86l and passive 
magnetic rods E2 can be found elsewhere. Here we focus on 
practical aspects of using DDM to characterise motile E. coli. 

We take a low-magnification movie of a cell population, and 
calculate the squared modulus of the Fourier transform (F.T., 
reciprocal vector q) of the difference of two images separated 
by r in time, and averaged over starting times: 

g(q,r) = (|F.T. [. T(x,y ; t + t) - I(x,y, t)f) . (7) 

If the image intensity I{x,y\ t ) at pixel position (x, y) at time t is 
linearly related to the cell density at the corresponding sample 
position and time, then g(q, r) is directly related to the density- 
density time correlation function, /(q, r), of the system, with 
q — |q| ~ 27r(irnage size) 1 to zrtpixel size) 1 being experimen¬ 
tally accessible. Indeed, under appropriate conditions, 

g{q, t) =A{q)\\ - f(q, r)] + B(q), (8) 

Sometimes called the normalised dynamic structure factor or 
intermediate scattering function (ISF), f(q, r) is in principle 
(but not in practice ||85l ) also measurably by dynamic light scat¬ 
tering. We assume isotropy, so that g and / only depend on q. 
In Eq.[8j A(q) is related to the form and (static) structure factors 
of the bacterial suspension, and B(q) is instrument noise. 

The ISF can be modelled once we know how cells move. For 
smooth-swimming E. coli in 3D mixed with non-motile cells, 

f(q, t) = (3e~ q I>T + ( \ - (j)e~ q Dt f P(v ) Sm ^ ^ dv, ( 9 ) 

J v qvr 


where P(v) is the speed distribution of the swimmers, and D 
is the free translational diffusion coefficient, which is mainly 
determined by the diffusivity of the non-motile cells (but see 
footnote[6]). We have ignored polydispersity in cell size (which 
induces a distribution in D) for simplicity. 

Self consistency demands that for straight swimmers mixed 
with diffusers, the dynamical quantities obtained by fitting Eq.|9l 
to data, in particular v = J () vP(v)dv, [3 and I), should be q- 
independent. Residual q dependence is nevertheless observed 
in practice. Fig. ib), e.g. due to tumbling] 1861: such q depen¬ 
dence is the main source of experimental uncertainties in the 
fitted parameters from DDM quoted in this work. 

The average body rotation angular frequency, O, can be de¬ 
termined using another high-throughput technique, dark-field 
flicker microscopy (DFM) ll89l . Essentially, the lowest peak in 
the Fourier transform of the time-dependent image intensity of a 
single cell in dark-field microscopy is O. In the low-Reynolds- 
number limit, v/Cl depends solely on the geometry of the swim¬ 
mers, and not on viscosity, provided that the medium is a New¬ 
tonian fluid (D- Thus, non-constancy of v /£1 indicates either 
changing cell geometry or non-Newtonian effects lf89l . 

Note that the relative merits of DDM and real-space tracking 
for obtaining v and [i has been discussed in depth before li86l . 

4.1. Imaging protocol 

We use d ~ 400 pm deep flat glass sample cells filled with 
* 150 pi of suspension and rendered air tight with petroleum 
jelly to prevent drift due to evaporation and stop replenishment 
of 0 2 after it is exhausted by respiring bacteria^] DDM shows 
that contact with petroleum jelly does not change (v,yS). 

Movies in phase-contrast illumination (PF10x/0.3 at 100 
frames/s, ~ 4000 images at 512x512 pixels => ~ 40 s movies) 
are recorded in an inverted microscope (Nikon Ti) with a 
Mikrotron high-speed camera (MC 1362) and frame grabber 


5 The OD depends on the wavelength used; we use 600nrn throughout. 
"Note that the bulk diffusion is enhanced by the swimmers (83J 


7 Calculations of the effects of tum bling on the ISF exist ED. 

8 See Sectionji] especially Fig.|lo[b), to see why we want no 0 2 supply. 
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Figure 4: Average swimming speed v (main) and non-motile fraction /) (inset) 
vs. time at three heights from bottom wall for WT AB1157 at 5 X 10 8 cells/ml. 


(Inspecta 5, 1 Gb memory) at 22 + 1°C. Custom LabVIEW 
recording software controlling the microscope stage allows 
scanning of many samples (giving good averages) repeatedly 
(giving time resolution), giving data such as Figs, pf] [9] and p~0] 

4.2. Observational complications 

The majority of experiments to date using synthetic swim¬ 
mers produced data at or near surfaces. Escherichia coli cells 
can swim in the bulk for extended time, although cells en¬ 
countering a wall tend to be trapped there (90l [911 . If per¬ 
formed sufficiently far from walls, DDM yields a genuine three 
dimensional speed distribution (85l [86l via fitting an ISF to 
Eq.0 Near a wall, the swimming is two-dimensional, and 
the sin (qvr)/qvT in Eq. [8] is replaced by J 0 (qvr), the zeroth- 
order Bessel function. However, it may not be possible to dis¬ 
tinguish between these functional forms unless P(v) is suffi¬ 
ciently narrow. Moreover, E. coli swims in circles (92), giving 
^/-dependent fitted dynamical quantities if the circle radii are 
comparable to the accessible range of 2ttlq. The complex mo¬ 
tion of cells ‘tethered’ to walls further complicate data analysis. 
We therefore typically image at 100 pm away from the bottom 
of our capillaries (cf. total cell length L < 10 pm). 

The sedimentation height of non-motile cells, Eq. [2] can be 
enhanced by 50% or more because swimmers enhance their 
diffusivity |83l . Nevertheless, l g ~ 0(10 pm), so that non¬ 
swimmers will accumulate on the bottom on a time scale of 
d/v s < 70 min (using v s » 0.1 pms" 1 , Eq. |4j. Measurement of 
WT swimmers 10 pm from the bottom capillary surface. Fig [4] 
inset, shows that /3 increases with time, until it saturates just 
over an hour (~ d/v s ) into the measurement. All non-motile 
cells have now sedimented into a layer of thickness ~ l g . Cor¬ 
respondingly, p is constant at any single bulk position (e.g., 50 
or 100 pm from the bottom, Fig.[4]inset) until the sedimentation 
front passes (at some fraction of d/v s ), whereupon p drops. 

Measurement of WT swimming E. coli in a sealed capillary 
at our standard height of 100 pm from the bottom shows that 
v decrease with time. Fig. [4] A priori, it is conceivable that 
this is because fast swimmers encounter surfaces and become 
trapped li90l[9l l more rapidly, thus leaving slower swimmers in 


the bulk. However, measurements at three heights, Fig. 0 give 
essentially identical results]/] so that the observed decreasing 
v(t) is not due to gradual kinematic accumulation at surfaces. 

Drift gives ballistic motions that masquerade as swimming. 
Evaporative drift can be eliminated by proper sealing, but un¬ 
avoidable drift occurs at short times due to loading. These tran¬ 
sients decay on a time scale > pL 2 p / for a liquid with viscosity 
j) and density p, where L is a characteristic length |93l . We used 
a 5 cmx8 mmx0.4 mm capillary throughout. The 8 mm dimen¬ 
sion leads to resolvable speed changes in our time window on 
the ss 1 min time scale, which probably accounts for the initial 
fast decay in v(t). Fig. [4] We find that decreasing this dimension 
gives faster initial decay. 

4.3. Fitting lore 

The decay of the ISF, Eq.[9] is due to diffusion (parameterised 
by D ) and ballistic motion (parameterised by v). Satisfactory 
fitting of the data in the form of g(q, r) using Eq.[9]and Eq.[8]to 
obtain D and v depends on convincingly decoupling these con¬ 
tributions, which is straightforward when the two decays have 
well-separated characteristic times. Fig. [3] Overlap of these two 
contributions occurs when their relaxation times approach, i.e. 
v' —* qD, which, for WT E. coli in bulk, D > 0.3 pm/s, and 
q < 2.2 pm -1 (using a lOx objective), occurs at v < 1 pm/s. 
However, due to the distribution of speeds, this limit is under¬ 
estimated. In practice, we can resolve v down to 2 to 3 pms -1 
in BMB, Fig. [5] Note also that if/I —> 0 or 1, the fitted value of 
I) or v» respectively will be subject to large errors. 

4.4. Motility plates 

When chemotactic strains are inoculated into ‘soft agar’ 
(> 0.2 wt%), they spread out in rings (94) . Non-chemotactic 
mutants m and the chemotactic WT when the agar is > 
0.3wt% f95) spread out instead as solid discs. The diameters 
of these rings or solid discs are often used to assay motility. 

We inoculated 5 pi drops of OD = 0.3 TB-grown E. coli into 
0.3wt% LB agar plates, measured the diameter of the solid 
colony discs after 22 h of incubation at 30 °C, and compared 
these with v measured from DDM in three different media. 
Fig. [5ja). For chemotactic WT cells, there is reasonably cor¬ 
relation between motility plates and DDM measurements in all 
three media, but especially glucose-supplemented BMB. 

The non-chemotactic spreading of growing cells can be mod¬ 
elled using the Fisher equation (%1 . This predicts a front speed 

of u ~ 2 V Do, where a is the growth rate, and D is the cells’ 

—2 

effective diffusivity. The latter is ~ v“r/3, with r being the 
average time between tumbles (971 . The colony diameter af¬ 
ter time At is d = uAt ~ 2vAt s[arF]3 oc v. Using the value 
of t ~ 1 s in bulk media (33l and the glucose-BMB data for 
HCB1, Fig. 0a), we find a ~ 0.0003 s 1 , or a doubling time of 
aT 1 In 2 > 30 min, which is reasonable under our conditions. 

Two caveats are in order. First, if the agar concentration al¬ 
lows chemotaxis, the dependence of the speed of the front (now 
a ‘ring’) on v becomes more complex |95j . Secondly, we find 


9 Note the longer time axis here than those shown later in Figs. 9 and 10 
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Figure 5: (a) Average swimming speed v (left axis) and (b) non-motile fraction 
/? for selected E. coli coli strains, from one batch culture in each case, before 
washing (in TB, white) and after washing (in BMB+1.5 mM glucose, grey or in 
BMB, black) at a cell density of 5 x 10 cells/ml. Red bars in (a) are diameters 
of solid, circular colonies in 22-hour LB motility plates (right axis). Note that 
for HCB437, the cells did not penetrate the agar but stayed on the surface. 


that HCB437, a smooth swimmer, did not penetrate the agar; 
its minimal spread on the plate surface does not reflect v in bulk 
media. Thus, motility plates must only be used with caution. 

5. Preparing motile E. coli 

We turn now to culturing bacterial swimmers. Readers need¬ 
ing to learn generic microbiological protocols should use stan¬ 
dard manuals l98l [99 ] and consult a microbiologist. We focus 
on extra procedures for obtaining motile cells. 

5.1. Culturing cells 

Our protocol is based on that of H. C. Berg (Harvard), and in¬ 
volves Luria broth (LB, originally ‘lysogeny broth ’ EES) and 
tryptone broth (TB). TB is a mixture of amino acids from hy¬ 
drolysing milk casein proteins; LB is tryptone broth plus yeast 
extract, the latter containing various carbohydrates not in TB. 

We grow single colonies from frozen stocks on LB agar 
plates at 30 °C overnight. A single colony is transferred from 
a plate to 10 ml of liquid LB and incubated overnight (~ 16 h) 
at 30 •C0 shaken (for aeration) at 200 rpm. Finally, cells are 
diluted 100-fold into 35 ml of TB and grown for ~ 4 h at 30 °C 
shaken at 200rpm. DDM shows that transferring to TB from 
LB results in higher v. As cell sizes (Table [T} and other proper¬ 
ties change with growth phase, it is important that growth times 
for different experiments should be kept constant. 


10 Growth is optimal at T = 37 °C, but lower T promotes swimming m\m. 
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Figure 6: (a) Average swimming speed v and (b) non-motile fraction /? during 
a washing process from TB to BMB (black) and from TB to TB (red) of WT 
AB1157, using filtration (filled) or centrifugation (open) at final cell density 
of 5 x 10 8 cells/ml. Inset: measured (symbols) and fitted (line) f(q, r) of these 
cells in TB before washing (turquoise, lower) and in BMB after 4 washing steps 
(purple, upper). Data sets are from one batch culture. 


5.2. Washing cells 

The next step is to transfer cells from TB to a ‘minimal 
medium’ with no exogenous nutrients. This prevents growth, 
which is an unwanted complication in the long ran. More im¬ 
portantly, we find that cells use oxygen very quickly in TB, and 
once oxygen is exhausted, v decreases dramatically. For WT 
AB1157 at OD = 0.3, this occurs within 10 min, which is not 
long enough for experiments. We therefore wash (by filtration) 
and transfer cells into Berg’s motility buffer (BMB), containing 
6.2 mM K 2 HP0 4 , 3.8 mM KH 2 P0 4 , 67 mM NaCl, and O.lmM 
undergo 4 successive filtrations (= 3 washes, see 
) until the concentration of TB is diluted 10 -fold 
by BMB. After the last wash, the suspension on top of the filter 
is removed and transferred to a 50ml centrifuge tube together 
with the filter, where cells are carefully suspended by rolling the 
tube. This yields 1-3 ml of highly concentrated cells (OD ~ 10), 
allowing the preparation of many dilute samples. 

Transfer from TB into BMB decreases v and increases /) for 
all strains studied. Fig. [5] This change in cellular motility is 
reflected in the ISF measured by DDM, Fig.[6ja) inset. In TB, 
the population of mostly motile cells (a single decay) has a nar¬ 
row speed distribution (oscillations). In BMB, two processes 
are clearly visible: a sizeable fraction of non-motile (diffusive) 
cells coexists with (ballistic) swimmers, whose speed distribu¬ 
tion is now considerably wider (no oscillations). 

Figure[6]also shows that washing by filtration has non-trivial 
effects on v and /), probably due to a combination of mechanical 
damage to flagella, regrowth (possible in TB) and change of 
medium. Pipetting may also have an effect. After 20 pipetting 
cycles in BMB using a 0.6 mm tip at normal expulsion speeds. 


EDTA. Cells 
Appendix D 
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Figure 7: Example of viable plate count vs. OD used for calibrating the latter. 

These values are typical for our culture protocol. 
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we find that v drops by ~ 50% and [3 increases by a factor of 
< 3. It is not until the tip diameter has been increased to « 
2 mm that the effect on (v,j8) begins to saturate. Such flow- 
driven damage of E. coli flagella has been noted before 1 1021 . 

If low cell concentrations are required or many strains need to 
be processed simultaneously, we use a bench-top centrifuge to 
wash small volumes (< 2 ml) 3 times at moderate speeds (2 min 
at 8000 rpm) to produce ~ 1 ml of OD « 0.3 cell suspension. 
Figure [6] shows that this procedure produces populations with 
higher (3 in TB, but it is considerably faster than filtration. 

5.3. Counting cells 

The most common protocol for relating OD to cell number 
is the ‘viable plate count’ 19811991 . Briefly, a sample is diluted 
to the relevant OD range (say, 0.05-1) using phosphate buffered 
saline (PBS)P~|One part of the sample is used to measure OD; 
another part is used to prepare tenfold serial dilutions (down to 
10 6 ) by transferring 100 pi sample to 900 pi PBS with vigor¬ 
ous mixing between each dilution. 100 pi of the final two dilu¬ 
tions are spread on replicate LB plates and incubated at 37 °C 
overnight. Next day, individual colonies can be seen. Each of 
these has arisen from a ‘colony forming unit’ (CFU), a short¬ 
hand for ‘probably but not certainly a single cell’. These are 
counted to arrive at the viable plate count in units of CFU/ml. 
The cell density of the original sample is then calculated taking 
into account the dilution. Plotting the OD vs the cell density 
from plate counting. Fig. [7] gives the desired calibration. 


6. E. coli energetics 


Self-propulsion requires energy. The energetics of synthetic 
swimmers is partially understood 1 1031 1041. Here we give a 
simplified introduction to E. coli bioenergetics. Details can be 
found in textbooks IH05lll06l[T07lll08lll09lfTT0ll (f] 

Bioenergy comes from transferring electrons from high- 
energy to low-energy bonds through a sequence of molecules. 


Figure 8: Schematic of E. coli bioenergetics, with high-energy compounds ATP 
and NADH picked out in bold. Glycolysis breaks down glucose (C 6 ) into two 
pyruvates (C 3 ), producing 4 ATP by substrate-level phosphorylation. Pyruvate 
is reduced to acetyl co-A with the emission of C0 2 . (Pyruvate is catabolised 
to acetate when glucose is abundant.) Acetyl co-A enters the tricarboxylic acid 
(TCA) cycle to produce 12 NADH per glucose. Each NADH donates it high- 
energy electron pair, 2 e _ , to respiratory enzymes in the inner membrane, which 
pass them to a terminal acceptor, here 0 2 , which is reduced to water, regen¬ 
erating NAD + (oxidised NADH) for glycolysis (dotted arrow). The drop in 
energy as 2 e _ passes along the respiratory enzyme chain is used to pump H + 
out of the cytoplasm, generating a proton motive force (PMF), which drives 
protons through rotary flagella motors and through F|F 0 ATPase to make ATP 
from ADP (chemiosmotic oxidative phosphorylation). 


The energy is used to pump protons out of the cell, giving rise to 
a proton motive force (PMF), which powers swimming. In aer¬ 
obic respiration, the electron donors are various reduced foods 
rich in C-H bonds, while the final electron acceptor is typically 
oxygen. The electron pair in the C-H bond lowers its energy as 
it is transferred (notionally as : H-) to an O-H bond (in water), 
where, because O is more electronegative than C, the electron 
pair interacts more strongly with positively charged nuclei than 
in C-H, and therefore has lower energy. The cell harnesses the 
energy released for chemical work. The ‘favourite’ supplier of 
high energy electrons for E. coli is or-D-glucose (C 6 H 12 Q 6 )P^| 
The full oxidation of this molecule in aerobic respiration gen¬ 
erates « 2.8 MJ mol -1 under physiological conditions 

OH 

+ 6 0, —■* 6C0 2 + 6H 2 0 AG » -10 3 k B T, 

( 10 ) 

Figure [8] summarises the aerobic catabolism (break down) 
of glucose, a C 6 compound. First, glycolysis via the Embden- 
Meyerhof-Parnas (EMP) pathway l4 | breaks glucose into two 
C 3 -halves, pyruvates (CH 3 COCOO ), a key metabolic interme¬ 
diate. Then, pyruvate is oxidised to acetyl coenzyme A (acetyl 



n PBS = 137 mM NaCl. 2. mM KC1, 10 mM Na,HP0 4 and 2 mM KH,P0 4 . 
12 An unusual, ‘coarse-grained’ introduction is given by Nobel laureate 
Christian de Duve fim . See also tutorials at fm) under ‘The Microbe’ tab. 


13 Glucose polymers constitute > 50% of dry terrestrial biomass H08l . 
14 When glucose is the sole carbon and energy source for E. coli, ~ 75% of 
it channels into the EMP pathway; the rest feeds into the hexose monophos¬ 
phate (HMP) pathway, which mainly provides precursors for biosynthesis, i.e. 
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co-A), another key metabolic intermediate, and fed into the ‘tri¬ 
carboxylic acid cycle’ (TCA = Krebs or citric acid cycle) to be 
further broken down. After glycolysis, C0 2 is produced with¬ 
out molecular 0 2 through the progressive enzymatic dehydro¬ 
genation of fragments of catabolised glucose. 

During glucose catabolism, the energy released is stored in 
high-energy compounds, mainly adenosine triphosphate (ATP) 
and (reduced) nicotinamide adenine dinucleotide (NADH). 
ATP is generated from adenosine diphosphate (ADP) by the ad¬ 
dition of an inorganic phosphate group: 

ADP + P ;—* ATP. (11) 

NADH comes from reducing NAD + , the oxidised form of this 
compound, partly via C0 2 -producing dehydrogenations: 

COOH C0 2 

NAD + + H + + 2 e“ ^—4 NADH. (12) 

The generation of 4ATP per glucose en route to pyruvate is 
known as ‘substrate level phosphorylation’ (SLP), where phos¬ 
phate groups are added to ADPs without using molecular 0 2 . 
The 12 NADH generated in glycolysis and the TCA cycle give 
rise to up to 96 ATPs through the process of ‘oxidative phospho¬ 
rylation’ . NADH donates its energetic electron pair to respira¬ 
tory enzymes in the inner membrane, ultimately passing them to 
0 2 as ‘terminal electron acceptor’, producing H 2 0 and regen¬ 
erating NAD + for glycolysis and the TCA cycle. Without such 
regeneration, e.g., because respiratory enzymes are poisoned or 
because 0 2 runs out, the TCA cycle falters. 

The energy released as 2 e drops down the potential energy 
ladder of respiratory enzymes pumps between 2 to 8 H + out of 
the cell (8 at high [0 2 ]) 11131 . generating a proton motive force 
(PMF) of 150 mV. Driven back to the cellular interior by the 
PMF, protons do work in active solute transport, turning flagella 
motor and making ATP in ‘rotary enzyme’ FjF 0 -ATP synthase 
(‘chemiosmosis’, requiring ~ 4 H ’ per ATP). 

When glucose is plentiful, it is incompletely oxidised to ac¬ 
etate (CH 3 COO ) via pyruvate: 


O 

// 

C 6 H 12 0 6 +20 2 - >2CH 3 C x +2C0 2 + 2H 2 0. (13) 

CTH + 

Such ‘overflow metabolism’ (or the ‘Crabtree effect’) generates 
less energy, and excretes acetate for potential later reabsorption 
and further oxidation when glucose is less plentiful 11.14 1 15 1. 

How E. coli adapts to decreasing [0 2 ] in its environment is 
complex 1 11611117 , 11 8 1 1. At [0 2 ] = 0 and without other ter¬ 
minal electron acceptors, E. coli ferments glucose into acetate, 
formate, lactate and succinate, using metabolic intermediates 
as electron acceptors and producing ATP by SLP, with acetate 


anabolism urn However, the HMP pathway can also feed pyruvate into the 
TCA cycle flQ9l . A fraction of TCA intermediates supplies anabolic precur¬ 
sors. Such ‘anabolic syphons’ together with overflow metabolism, reaction [TT] 
mean that the stoichiometry in reaction! lOlis never satisfied on average. 
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Figure 9: (a) Average swimming speed v vs. time for the capillary protocol in 
BMB at cell densities indicated in the legend (cells/ml). Lines guide the eye. 
The ‘crash’ in v in each case is due to 0 2 exhaustion. Data sets from one batch 
culture, (b) Corresponding 0 2 exhaustion time as a function of inverse of cell 
density. Error bars relate to the accuracy of measuring the ‘crash’ time. 


production generating most ATP 11161 119 1. In such ‘mixed- 
acid fermentation’, respiratory enzymes are inactive, and other 
pathways maintain the PMF, including F, F (l -ATPase working 
in reverse, hydrolysing ATP to ADP and using the energy to 
pump protons out of the cell 1107 IK) 1171 . The PMF powers 
flagella motors. Within a wide range, the motor speed is pro¬ 
portional to the PMF 11201 . If the PMF is short circuited (e.g. 
by making the membrane permeable to protons) cm or drops 
to zero when ‘fuel’ runs out, swimming ceases. 

Under starvation conditions (i.e., no exogenous nutrients), 
such as occurs in BMB, E. coli obtains energy by metabolis¬ 
ing intracellular resources. Such ‘endogenous metabolism’ is 
complex li 122l 1123 ,1 and poorly understood. For starved cells 
grown in TB, catabolism of free amino acids and free or RNA- 
derived ribose (a sugar) feeds into the TCA cycle 11081 . so that 
0 2 or other terminal acceptors are needed. 


7. E. coli swimming powered by endogenous metabolism 

Escherichia coli and other bacteria 1124: , 12511126 1 can swim 
powered entirely by endogenous metabolism. We now probe 
the time dependence of v. All data are obtained from DDM with 
WT AB1157 in BMB prepared using our standard protocol. 

7.1. Observations and practical implications 

Figure [9j a) shows the average swimming speed as a function 
of time, v(t), in sealed capillaries at various cell densities, n. 
Cells slow down, approximately linearly with time. If ~ 10% 
slowing down is tolerable, then there is a useable window of 
« 1 h at n — 2 x 10 9 cells/mL, which ends in a ‘crash’. 

The time of this ‘crash’, t c , drops with n\ we suggest that it 
is due to 0 2 depletion. Quantitatively, nt c ~ constant, Fig.|9jb), 
as expected if the initial [ 0 2 ] and the 0 2 consumption rate per 
cell, Q, are invariant with n. Assuming that the BMB was ini¬ 
tially saturated with 0 2 , and using literature solubilities fll27l . 
which are only very weakly dependent on ionicity II 1281 . we 
estimate Q at each cell density. Fig. |9]b) inset, which is 
weakly dependent on n, rising from Q endo ~ 2 amol/min/cell 
at n = 1.67 x 10 9 cells mf 1 to Q endo < 4 amol/min/cell at 
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Figure 10: Effect of glucose on time-dependent average swimming speed v(t) 
for WT AB1157 at 5 x 10 8 cells/ml. (a) At low [Glu], glucose depletes first, 
causing the rapid drop (III) in v. Inset: a longer run at [Glu]=0.07 mM reveals a 
subsequent sharp drop (IV) due to 0 9 depletion, (b) At high [Glu], 0 2 depletes 
first, causing the first sharp drop (V) in v. The second, larger sharp drop (VI) in 
v is due to glucose depletion. Datasets are from one batch culture, (c) Deple¬ 
tion time of glucose in aerobic condition versus glucose concentration. Inset: 
aerobic glucose consumption rate versus [Glu]. (d) Speed v versus molecule 
concentration normalised to the speed in BMB, v 0 , for three small molecules, 
measured « 5 min after mixing and filling the capillary. Filled symbols are con¬ 
centration markers for Fig.[l2] Error bars in (c) and (d) are standard deviations 
from averaging data from two and three independent batch cultures. 


n = 5.55 x 10 9 cells ml -1 (‘a’ = ‘atto’ = 1CT 18 ). These val¬ 
ues are comparable to those measured before using a different 
strain and under somewhat different conditions B1221lp^] 

The highest density data in Figs. [9] were obtained at n « 
5 x 10 9 cells/ml, or <f> ~ 0.5% of cell bodies. Extrapolating from 
Q at this n, t c « 1 min at </> « 5%; this is an upper bound, since 
Q probably increases with n. Oxygen exhaustion therefore pre¬ 
cludes experimentation at </> > 1% under these conditions. 

7.2. Physiological interpretation 

The observed slowing down in v at t < 1 h is probably not 
due to substrate exhaustion, as extrapolated literature data [ 1221 
1123 1 suggests this happens in ~ 10 h. That dv/dt is identical at 
different n. Fig. [9ja), supports this suggestion. 

The ubiquitous signalling molecule cyclic di-guanylate (c-di- 
GMP) can slow down bacterial rotary motors by binding to the 
protein YcgR to form a molecular ‘brake’ H129II . However, we 
found that a A YcgR mutant behaved identically to the WT in 
sealed capillaries; so c-di-GMP is unlikely implicated. Instead, 
we infer a slow decrease in the PMF driving the rotary motors. 
The adaptive value and mechanism of this putative starvation 
response is unknown, but it can affect many PMF-dependent 
cellular functions, e.g. ATP synthesis by F^Q-ATPase. 


15 &ndo ~ 10-20 ml/h/mg dry cell mass (m -> 2 -4 amol/min/cell using 
0 2 molar volume = 24 dm 3 and dry cell weight « 0.3 x 10 -9 mg/cell (Ml. 



Figure 11: Simultaneous measurements of average swimming speed v (o, left 
axis) and [0 2 ] (□, right axis), using fluorescence lifetime imaging microscopy 
(FLIM) vs. time in a sealed capillary for WT AB1157 at « 10 9 cells/ml and 
[Glu]=0.5 mM. For FLIM, 47 mM of the 0 2 -sensitive dye ruthenium tris(2,2’- 
dipyridyl) dichloride hydrate (RTDP) was added to BMB (which did not affect 
the cells’ motility). RTDP was excited using sub-ps pulses (T = 450nm, ~lmW 
at 1MHz repetition rate) and its fluorescence imaged using a gated intensified 
CCD camera (Picostar HR-12QE, La Vision GmbH, Germany). The data was 
fitted with a single exponential decay yielding a (homogeneous) fluorescence 
lifetime map, which was averaged to estimate [0 2 ]. Inset: v at 5x 10 8 cells/ml to 
better show the saturation before 0 2 depletion. Data sets are from independent 
batch cultures with v{t = 0) « 12 pm s -1 and « 15 pm s -1 , respectively. 


8. E. coli swimming powered by glucose 

If glucose is present, E. coli utilises it first, and suppresses 
the expression of enzymes for processing other nutrients 11301 . 

Glucose increases the swimming speed of E. coli nnn (see 
also Fig.[5ja)) and other enteric bacteria Il3ll . We now demon¬ 
strate how to use glucose to enable E. coli AB1157 (cultured by 
standard protocol) to swim at a constant v for extended time. 


8.1. Observations and practical implications 
Figure 
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a) and (b) show v at n = 5x10 cells/ml of WT 
AB1157 sealed in a capillary in BMB supplemented with low 
and high glucose concentrations respectively. At 0.03 mM < 
[Glu] < 0.07 mM, Fig.| 1 0[a), v rises gradually with time (II) af¬ 
ter initial transients (I), until it drops suddenly (III). Thereafter, 
it decreases slowly with time. A longer run at [Glu] = 0.07 mM, 
Fig. pi)]; a) inset, reveals a second rapid drop (IV) that is sharper 
and of larger amplitude. The long, slow decrease with time and 
the final sudden drop are reminiscent of ‘endogenous propul¬ 
sion’, Fig. [9] Thus, the data suggest that glucose is consumed 
aerobically until it is depleted, whereupon cells revert to en¬ 
dogenous metabolism until O, is, in turn, depleted. 

Consistent with this, the time of the first sharp drop (III) is 
proportional to [Glu], Fig. [TT)|c). and the deduced specific glu¬ 
cose consumption rate, inset, is nearly independent of [Glu], 
averaging to G ~ 6.9 ± 0.5 amol/min/cell. We also measured G 
over a range of [Glu] biochemically (Appendix E i, and found 
a comparable value of G « 5 ± 0.6 amol/min/cell. 

Figure |T0[c) implies that there will be a [Glu] at which 0 2 
depletion will occur before glucose depletion. This indeed hap¬ 
pens at the higher [Glu] explored in Fig. [T0}b). Now, the very 
sharp, first ‘crash’ (V) in v is due to O, depletion, at which 
point there is still plentiful glucose in the medium. Consistent 
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Figure 12: Normalised average swimming speed, v/v 0 , from Fig 1 1 0| b) as a 
function of the osmolarity. Note the bottom split axis with a linear scale from 
150 to 190 mOsm (bottom left) and a log-scale from 200 to 2000 mOsm (bottom 
right). Filled symbols corresponds to markers in Fig.|10[d). 


Figure 13: Average steady-state swimming speed v of E. coli MG1655 (n = 
10 s cells/ml) expressing SAR86 y-proteorhodopsin at different intensities of in¬ 
cident green light (,t = 510 - 540 nm) after 0 2 depletion; data obtained starting 
from the highest intensity. 


with this interpretation, the time of the first abrupt speed de¬ 
crease (V) is independent of [Glu], Again using literature 0 2 
solubilities 1 1271 , 1281. we find an average 0 2 utilisation rate of 
Qciu = 22.6 + 0.6 amol/min/cell, which, as expected, is consid¬ 
erable higher than the Qendo ~ 2-4 amol/min/cell reported by us 
(Section |7TT| i and others 11221 for endogenous metabolism. 

Independent evidence for 0 2 depletion comes from monitor¬ 
ing the fluorescence life time r of an 0 2 -sensitive dye Iil32l . 
By measuring r and using our experimentally derived calibra¬ 
tion [0 2 ] = 568.1 |iM (604ns/r - 1) (consistent with 111331 ). we 
can assay [0 2 ] in our sealed samples. Fig. 11 which drops to 
~ 0 precisely at the point of the abrupt decrease (V) in v. 

Subsequent swimming powered by anaerobic glucose 
metabolism shows v constant to ±1 pm s 1 until glucose deple¬ 
tion and v abruptly drops by a factor of 3 or more. Fig. 1 10(c) 
(VI). Thus, cells in initially aerated BMB supplemented with 
enough glucose will swim at approximately constant average 
speed after 0 2 depletion and before glucose depletion. 

The period of constant-v swimming can, at first sight, be ar¬ 
bitrarily extended by increasing [Glu], Fig. [TOjb), up to sat¬ 
uration (~ 5M at 25 °C). However, a practical limit exists. 
Figure p~0|d ) plots w 0 , the average swimming speed at vari¬ 
ous [Glu] normalised to the average speed at [Glu] = 0, both 
measured immediately after loading into sealed capillaries. At 
[Glu] < 10 5 mM, glucose depletion occurs before measure¬ 
ments could begin, so that w 0 = 1. At [Glu] > 50 mM, w 0 
catastrophically drops. Repeating using L-lactate, which E. coli 
can metabolise, and sucrose, which it cannot, also shows the 
same abrupt drop, which we believe is due to osmotic shock. 

Figure [TJJreplots our data against solution osmolarity (deter¬ 
mined using an Osmomat30, Genotec, Germany), showing that 
v always starts to decrease at an osmolarity between 0.2 and 
0.4 Osmom, comparable to the osmolarities needed to start de¬ 
creasing cell volumes in previous studies 1 1341 . Thus, there is a 
limit to the highest useable [Glu] for maintaining constant vp*] 


16 As expected from previous work fl34l . the results in Fig. EH were time 
dependent: waiting longer produced partial recovery in the ‘crashed’ speed. 
Nevertheless, a definite upper limit exists at« 100 mM small-molecule solutes. 


8.2. Physiological interpretation 

Reaction [T0]for complete glucose oxidation requires an oxy¬ 
gen to glucose stoichiometry of Qq\ u /G = 6. We measure 
(2 gi u /G ~ 3.3 + 0.33, i.e. about 50% of the glucose consumed 
is fully oxidised. A combination of overflow metabolism and 
‘anabolic siphon’ (footnote p~4]) probably explains this finding. 

Figure m shows the early-time portion of v(t) on an ex¬ 
panded time axis. After initial transients (I), v(t) rises (II) and 
then saturates (this is particularly clear in the inset) before the 
crash (V) due to 0 2 exhaustion. Immediately before these ex¬ 
periments, our cells have been harvested from TB, which con¬ 
tains little carbohydrates. While the operon of genes neces¬ 
sary for glucose transport ( pts ) is constitutively expressed irre¬ 
spective of whether glucose is present, exposure to glucose can 
increase the expression of some pts genes by up to threefold 
1 135l[ 136 1, which probably causes the rise (II) of v(t). 

Averaging over data sets, the ratio of the average speeds be¬ 
fore and after the 0 2 ‘crash’ (V) in the range 1.3 to 1.6. The 
PMF of K-12 E. coli (strain AN387) grown aerobically on glyc¬ 
erol and fermentatively on glucose was previously measured to 
be -160mV and -117mV respectively 1113711 : 160/117 = 1.36 
is remarkably close to our speed ratio. 

8.3. Elow much energy does swimming use? 

The detailed mechanics of a peritrichous swimmer (= many 
flagella distributed over the cell body) is not yet understood 
E49I , so that the vast majority of the literature on E. coli is 
couched in terms of propulsion by a single ‘effective flagellum’. 
The power developed by the ‘effective motor’ turning this ‘ef¬ 
fective flagellum’ has been estimated to be P « 4 x 10~ 16 W for 
swimming in minimal medium 11381 . Each H + driven into the 
cell by a PMF of « 150 mV U37| can perform w = 2.4 x 1(T 20 J 
of work. If we assume that flagella motors extract work with 
100% efficiency from the protons, then P/w < 2 x 10 4 H 1 s 1 
is needed to power the ‘effective motor’. In reality, of course, 
a somewhat higher flux would be needed, because the motor is 
not 100% efficient in harnessing energy from the PMF. 

On the other hand, full oxidation of one glucose gives 12 
NADH; each NADH can export up to 8 H + , Thus, G/2 ~ 
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Figure 14: A non-raotile Keio knockout ED, A(fLiC)::kcin, where the flagella 
protein gen efliC (red) is replaced by the kanamycin resistance gene kan (blue), 
flanked by DNA sequences (purple) that allow its easy excision. 


4 amol/min/cell = 4 x 10 4 /s/cell fully-oxidised glucose 

molecules export up to 384 x 10 4 H + /s/cell. This is in large ex¬ 
cess of the ~ 10 4 H + /s/cell we estimate as needed to power the 
~ 50% increase in v upon addition of glucose. Thus, swimming 
accounts for only a small part of a cell’s energy budget. 

8.4. Other small molecules 

We have also studied swimming powered by acetate (C 2 ), 
glycerol and lactate (C 3 ), xylose (C 5 ), galactose (C 6 ), and mal¬ 
tose and lactose (C 12 ), which can be expected to generate dif¬ 
ferent amounts of energy per molecule. Interestingly, however, 
the maximum increase in speed in each case is approximately 
the same as observed for glucose (about 50%, see Fig. [To} d) 
for L-lactate). This suggests that the flagellar motor speed satu¬ 
rates at high PMF (as observed for Bacillus subtilis If 139 l). and 
supports our conclusion that swimming uses only a small frac¬ 
tion of metabolic energy. At high enough concentration of these 
alternative carbon sources, we find that the post-0 2 -exhaustion 
speed drops to zero (for lactate, cf. 11401 ), rather than to some 
intermediate, constant level as observed for glucose. This is 
because E. coli cannot ferment these substrates. 


9. E. coli swimming powered by light 


Some micro-organisms directly harvest energy from light. 
Thus, photosynthesizers use photons to split water to generate 
high energy electrons. Alternatively, light is harvested directly 
to pump protons using proton pumps such as proteorhodopsin 
(PR). PR, first discovered through the genomic sequencing of 
marine bacteria 11421 . has been incorporated into E. coli 11431 . 

Most PRs identified to date can generate PMF up to or 
slightly above WT levels mm To enable illuminated PR to 
take over PMF generation, other sources need to be deactivated. 
This has been demonstrated on single cells using poisons for the 
respiratory enzymes (cf. Fig. [8]) or by 0 2 exclusion. As the mo¬ 
tor, and therefore flagellum rotation speed, is proportional to 
the PMF, illuminating poisoned cells expressing PR leads to a 
controlled increase in flagellum rotation speed 11431 , 

Here, we demonstrate this effect in a cell population using 
MG 1655 transformed with plasmid pBAD-HisC-PR express¬ 
ing PR when induced by arabinose (DM-4). We already know 
that 0 2 is exhausted after a period when swimmers are sealed in 
capillaries. Fig. [9] After our transformed cells have exhausted 
0 2 , their swimming is wholly powered by PR-generated PMF. 
The initial increase in speed with incident intensity. Fig. 13 sat¬ 


urates at high intensities, in agreement with previous measure¬ 
ments of the rotation rate of individual flagellar motors 1 1431 . 
This and other PR mutants therefore enable simple external 
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(kanamycin resistance) 
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Figure 15: Schematic of transposon mutagenesis. See text for details. 


control of their swimming speed, making it feasible to change 
the samples activity both in time and space. They are living 
analogues of light-driven synthetic swimmers 191 11451 . 


10. E. coli - the genetic toolbox 


A great attraction of E. coli is the availability of many mu¬ 
tants and a versatile genetic ‘tool box’. Thus, e.g., the FliC 
flagella protein could be mutated to include a cysteine residue to 
make disulphide bonds with dyes, or the cheY chemotaxis gene 
could be deleted to give smooth (non-tumbling) swimmers. We 
briefly introduce two ‘off the shelf’ mutant libraries and two 
methods for mutagenesis to illustrate principles and introduce 
terminology. Details can be found in the textbooks, e.g. ons. 

A major resource is the Keio collection 11411 of ‘sin¬ 
gle knockout’ mutants of strain BW25113 (closely related to 
MG 1655 but more amenable to molecular biology), available 
on a non-profit making basis to academics. Each non-essential 
gene was precisely deleted and replaced with an excisable 


kanamycin resistance gene. Fig. 14 


Many mutagenetic methods use plasmids, extrachromoso- 
mal pieces of DNA that can be ‘transformed’ into bacteria and 
maintained independently. Figure [l5]illustrates one strategy to 
carry out ‘transposon mutagenesis’. It relies on a temperature- 
sensitive plasmid, which is only maintained within the cell at 
the ‘permissive’ temperature, and is lost at a (higher) ‘restric¬ 
tive’ temperature. This ‘suicide’ plasmid acts as a ‘vector’ (car¬ 
rier) for a transposon, a piece of genetic code that can ‘jump’ 
between DNA locations. Here, the plasmid carries transposon 
Tn5, which confers resistance to the antibiotic kanamycin. 

When the transposon jumps, or transposes, from the plasmid 
onto the bacterial chromosome, it may insert into a gene, gener¬ 
ating a mutation. Following a shift of the E. coli from the per¬ 
missive to the restrictive temperature, the plasmid will be lost. 
Then, all kanamycin-resistant cells carry a chromosomal copy 
of the Tn5 element transposed onto the chromosome, where it 
can be stably maintained. Crucially, the transposon jumps into 
random chromosomal locations, and so can generate mutations 
in any gene. A particular mutation can be selected by a change 
in phenotype and the insertion site mapped by sequencing. 

PI transduction is another popular mutagenesis technique. 
Here DNA is transferred from a donor to a recipient cell us- 
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Figure 16: Schematic of PI phage transduction. (T) Phage attaches to cell. © 
Phage injects DNA (black) into cell and uses its machinery to replicate as well 
as fragment cell’s genome (red). @ DNA is packaged into viral particles, some 
of which contain fragments of the cell’s DNA, giving a ‘Trojan’ phage. (4) Cell 
lyses, releasing all viral particles. © ‘Trojan’ phage infects target cell. © 
Donor DNA fragment is incorporated into target cell genome by recombination. 


ing the PI bacterial vims (a bacteriophage). When PI infects 
a host cell, it not only packages its own genome into the vims 
particle, but also can package fragments of the host chromo¬ 
some. When these vims particles are purified and used to infect 
the target strain of E. coli, DNA homology between the donor 
DNA packaged in the virus particle and the recipient chromo¬ 
some can generate insertion of the donor DNA. 

More generally, plasmids constitute a vital resource for in¬ 
troducing new genes into E. coli , e.g. transforming with a plas¬ 
mid carrying the GFP gene will render cells fluorescent. Each 
member in another useful library, the ASKA collection 11471 . 
consists of an E. coli strain containing a plasmid in which a 
single gene from E. coli strain W3110 has been cloned. The 
DNA encoding each gene has been fused to sequences that per¬ 
mit precise temporal control of the expression of the associated 
protein and its fluorescent tagging with GFP. 

We end with two ‘health warnings’. First, genetic modifica¬ 
tion in motility-related work typically starts with parent strains 
selected for motility (e.g. HCBxx or RPxx, Table |B.2) . Many of 
these retain wild-type features not optimised for DNA manip¬ 
ulation. Moreover, the genome of many ‘motility favourites’ 
remains less well characterised compared to ‘canonical’ strains 
such as MG 1655 |26l . Sequencing and annotating the genome 
of some of these strains is a desirable future step. 

Secondly, care must be taken in obtaining a good reliable 
starter culture of the strain that one chooses to work with, ei¬ 
ther from the originating author or from a national culture col¬ 
lection, e.g. ATCC mi and CGSC EH in America; NCIMB 
1150) in the UK; DSMZ Q7Q in Germany; and CGMCC in 
China 115211 . Bacterial genomes are changeable, so that strains 
deposited at the back of freezers twenty years ago whose prove¬ 
nance is uncertain are not good starting points. 

11. Summary and conclusions 

Escherichia coli is increasingly used in active colloids ex¬ 
periments as an alternative to or to contrast with synthetic col¬ 
loidal swimmers. We have reviewed, explained and proposed 
a set of procedures that, together, should enable the collection 


of data from such experiments between different laboratories 
that can be quantitatively compared. At the heart of these pro¬ 
cedures is the use of a new, high-throughput method to charac¬ 
terise E. coli motility, differential dynamic microscopy (DDM), 
which delivers the average swimming speed and motile fraction 
averaged over ~ 10 4 cells in a matter of minutes 11851 1861. Us¬ 
ing DDM, we have established that E. coli swims at constant 
speed when motility is powered by mixed-acid fermentation, 
Fig. |T0[b). Along the way, we have explored the use of DDM 
to interrogate aspects of cellular physiology via the dependence 
of swimming speed on PMF. Our results suggest that E. coli 
cells powered by endogenous metabolism actively reduce their 
PMF with time. It appears that about 50% of the glucose un¬ 
der our conditions is used for generating energy, and that a very 
small fraction (maybe ~ 1%) of this energy is used for motility. 

Whether our proposed procedures, or something similar, are 
followed, it is important that sufficient information is recorded 
in publications using E. coli as a model active colloid to enable 
faithful reproduction of experiments. We therefore propose a 
‘checklist’ of what needs to be reported in | Appendix A 

More needs to be done before E. coli can become an ideal 
model for active matter experiments. Biologically, it would be 
desirable to sequence the genomes of some favourite motility 
strains. Armed with accurate genomic information, it should be 
possible to design bespoke swimmers that have ideal character¬ 
istics for different aspects of active matter research, e.g., a strain 
that does not depend on either poison, or 0 2 depletion, or light 
to display constant swimming speed over extended time. Theo¬ 
retically, it would be useful to construct coarse-grained models 
that take into account the flagella bundle. Predictions from such 
theoretical models confronted with data from experiments us¬ 
ing increasingly well-characterised E. coli suspensions should 
propel the subject of active colloids significantly forward. 
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Appendix A. Protocol reporting checklist 

To facilitate reproducing and comparing results, we suggest 
that published methods should always include: 

• Strain name and details of any genetic modifications 

• Growth conditions (temperature, growth medium compo¬ 
sition, cell density of harvest) 

• Washing procedure (method: filtration or centrifugation, 
number of washing steps, exact composition of motility 
medium) 

• Conditions for motility experiments (cell concentration 
and how determined, sample cell details: material and 
dimensions, sealing, position and conditions for im¬ 
age/video acquisition) 

• Details for image analysis and data processing 

Appendix B. Strains 

Table lB.2l summarises a number of strains of E. coli that have 
been used by us and others in active colloids and related motil¬ 
ity work. Note that publication in a journal carries the obli¬ 
gation of curating and sharing the strains used in a published 
work.However, one should not give away a strain obtained from 
another laboratory; instead the protocol is to direct enquiries to 
the strain originator. Note that intra-strain variability can arise 
rapidly because certain genetic elements (called insertion se¬ 
quences, IS) can ‘jump’ into regions of the genome controlling 
flagella gene expression, enhancing the motility [? ]. 


Appendix C. Nomenclature for bacterial genetics 


Researchers working with bacteria, especially if they collab¬ 
orate with biologists, need some facility in the shorthand used 
to described genotypes [ |157 , f58l , which in turn requires some 
knowledge of molecular genetics Ea. Here we give a very 


brief introduction based on ‘worked examples’ from Table B.2 


Each gene is given a three-letter, lower-case italicised name 
(possibly followed by an extra capital letter for different mem¬ 
bers of a gene family), often suggestive of its function; differ¬ 
ent versions of a gene (alleles) are distinguished by numbers. 
Thus, tsrl 02 1 is allele 7021 of the tar gene in E. coli coding 
for part of a receptor system for sensing L-serine and related 
amino acids. Deletion of this allele in HCB427 is indicated 
by A(tarp()2\. The protein product of a gene is denoted by 
the same three letters unitalicised and capitalised. Thus, in the 
fliC( S353C) mutation carried by ADI, residue 353 of the FliC 
protein used to build flagella filaments has been mutated from 
serine (S) to cysteine (C)p] Sometimes, a different gene is in¬ 
serted in place of a deleted gene. Thus, in JSL1, a cassette 
of genes conferring resistance to the antibiotic kanamycin has 
been inserted in place of the deleted cheY gene for chemotaxis 


17 See the textbooks, e.g. 11061 . for single-letter abbreviations of amino acids. 


- A(cheY::kan); the advantage is that such mutants can be eas¬ 
ily selected by growth in the presence of kanamycin. Finally, 
many bacteria carry plasmids, such as plasmid pHC60, which 
carries tetracycline resistance (Tet R ) and expresses green fluo¬ 
rescence protein (GFP) constitutively (i.e. continuously rather 
than subject to environmental control). 

Genotype of strains can be found from various online 
databases, including the various culture collections 1 1148 . 1491 
11501115 T , 1521 . Information on individual genes are available 
from appropriate databases, e.g. cm 


Appendix D. Filtration 

To transfer the bacteria from TB growth media to BMB, we 
use a sterile Nalgene filtration unit consisting of two compart¬ 
ments separated by a Milipore 45 pm HATF filter II1 601 . The fil¬ 
ter is soaked in BMB for ~ 5 min before being placed centrally 
on the unit with sterile tweezers. After screwing tight the up¬ 
per and lower halves of the unit, 35 ml of cell culture is slowly 
poured onto the filter. A tap-powered suction pump is used to 
enhance the flow rate, with a 70% Ethanol bath between the 
unit and the suction tap to prevent water supply contamination. 

When the level of TB had fallen to w 3 mm above the filter, 
leaving » 3 ml of liquid, 35 ml of BMB is pipetted into the unit. 
Ensuring that the filter does not run dry during the procedure 
minimises the number of non-motile cells. This washing step is 
repeated 3 times. The water flow is adjusted to maintain a con¬ 
stant filtration rate so that each wash step takes 10 min. When 
~ 1 to 3 ml of filtrate remains after the third wash, a sterile 
cut pipette tip is used to transfer the filtrate to a plastic 50ml 
polystyrene test tube (Greiner). (Smaller tubes are not wide 
enough for next step.) The filter is removed from the unit using 
sterile tweezers and deposited on the side of the tube. By gen¬ 
tly rolling the liquid over the filter, bacteria were resuspended 
to reach a final OD between 5 and 15, depending on the volume 
of liquid left on-top of the filter after the final washing step. 

Appendix E. Glucose assay 

We used a Sigma GAGO20 glucose assay kit. Hydrogen per¬ 
oxide from the oxidation of glucose by glucose oxidase reacts 
with o-dianisidine in the presence of peroxidase to form a col¬ 
ored intermediate, which further reacts with sulfuric acid to 
form a more stable, pink final product. The OD measured at 
540 nm is proportional to the original glucose concentration. 
Calibration runs showed that in our equipment, the linear range 
for OD 540 vs glucose concentration extended to ~ 600 pM. 

Bacterial cells were grown using the standard motility pro¬ 
tocol, washed and OD adjusted to 0.3 as for motility exper¬ 
iments, except that the buffer here lacks EDTA, which inter¬ 
fered with glucose oxidase. Cells and glucose were mixed and 
samples taken at different time points, put on ice and filtered 
(0.22 pm pore size) as quickly as possible. The filtered super¬ 
natant was used according to manufacturer’s protocol to mea¬ 
sure glucose concentration. Glucose utilization rates were cal¬ 
culated using glucose and time difference as well as wet or dry 
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Table E.3: Glucose utilisation rate 


Initial [Glu] (mM) 

rate ((xmol/min/g wet cells) 

0.50 

4.9+ 1.0 

0.20 

5.2 ± 0.9 

0.10 

6.1 ± 1.6 

0.05 

4.7 ± 0.8 


weight of cells in 1ml of OD = 0.3 suspension. Results are re¬ 
ported as averages and standard deviations using different time 
periods measured in one set of experiments. There is no sys¬ 
tematic trend. We therefore average the value determined at 
different [Glu] to give an estimated glucose utilisation rate of 
5.2 ± 0.6 pmol/min/g wet cells, or ~ 5 ± 0.6 amol/min/cell. 
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Symbol 

Typical units 

Brief definition 

[■•J 

pM or mM 

concentration of... 

x or ( x ) 

- 

average of x 

Latin 



c 

pM or mM 

small molecule cone. 

D 

[xm 2 s -1 

diffusion coefficient 

d 

mm 

colony diameter 

G 

amol/min/cell 

glucose consumption rate 

8 

rns^ 2 

gravitational acceleration 

L 

[xm 

cell body + flagellum length 

1 

[xm 

cell body length 


[xm 

gravitational length 

n 

cells/ml 

cell density 

n 

- 

refractive index 

Q 

amol/min/cell 

0 2 consumption rate 

q 

pm -1 

scattering vector 

T 

Kor °C 

temperature 

t 

s or min 

time 

*c 

s or min 

time of speed crashes 

u(t) 

- 

m/v 0 

V 

pm J 

volume 

V 

pm s" 1 

average speed 


pm s 1 

o 

IK 

III 

o 

II 

v s 

pm s” 1 

sedimentation speed 

W 0 (c) 

- 

w ; o(c) = v 0 (c)/v 0 (c = 0) 

Greek 



a 

s-' or min -1 

growth rate 

p 

- 

non-motile fraction 

y 

- 

l/cr 

n 

mPas 

viscosity 

A 

jim 

wavelength in vacuo 

p 

pmcmV -1 s i 

electrophoretic motility 

V 

- 

number of protons (H + ) 


pN s pnT 1 

friction coefficient 

p 

gcnT J 

(mass) density 

cr 

{xm 

cell body diameter 

T 

s 

delay time, run time, life time 

$ 

- 

cell body volume fraction 
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Table B.2: E. coli strains and plasmids 


Name 

Motility genotype/relevant characteristics 

Source 

E. coli strains 



MG 1655 

motility WT 

Laboratory strain 1261 

AB1157 

motility WT 

Laboratory strain [ 1531 

RP437 

motility WT 

J. Parkinson 11541 

HCB1 

motility WT 

H. C. Berg |155| 

BW25113 

Parent strain of Keio collection, motile 

1141] 

HCB437 

A(tsr)102\ A(trg)l 00 A(cheA-cheZ )2209 

Defective in chemotaxis, a smooth swimmer 

H. C. Berg 

DM4 

MG1655 + proteorhodopsin-expressing pBAD-HisC-PR 

This work 

JSL1 

AB1157 A(cheY::kan) 

Smooth swimmer 

This work 

ADI 

AB1157^/C(S353C) 

Motility WT with FliC mutation to bind to Alexa dye 

This work 

AD2 

ADI A(cheY::kan) 

Smooth swimmer with FliC mutation to bind to Alexa dye 

This work 

AD 3 

ADI A(ycgR::kan) 

Motility WT, unable to ‘brake’ 

This work 

Plasmids 



pHC60 

Tet K , GFP constitutive 

1156] 

pBAD-HisC-PR 

Amp K , arabinose inducible proteorhodopsin 

034] 
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